Isotropic Cosmological Singularities: 
Other matter models 



K.P.Tod 

Mathematical Institute 
Oxford 

February 7, 2008 



Abstract 

Isotropic cosmological singularities are singularities which can be 
removed by rescaling the metric. In some cases already studied ([H , 0) 
existence and uniqueness of cosmological models with data at the sin- 
gularity has been established. These were cosmologies with, as source, 
either perfect fluids with linear equations of state or massless, collision- 
less particles. In this article we consider how to extend these results 
to a variety of other matter models. These are scalar fields, massive 
coUisionless matter, the Yang-Mills plasma of j|] and matter satisfying 
the Einstein-Boltzmann equation. 



1 Introduction 

An isotropic cosmological singularity is one which can be removed by rescal- 
ing the spacetime metric with a single function which becomes singular on 
a smooth spacelike hypersurface, say S, in the rescaled space-time 
In the rescaled but unphysical space-time the curvature is finite at S. The 
Weyl curvature with indices suitably arranged is conformally- invariant, so 
that it is also finite at the curvature singularity in the unrescaled physical 
space-time. The motivation for studying this class of singularities comes 
from Penrose's Weyl Curvature Hypothesis (|^], |^]) which can be inter- 
preted as the hypothesis that the Weyl curvature must be finite at any 
initial singularity. 

In order to make mathematical progress with the Weyl Curvature Hy- 
pothesis one must first say exactly what is meant by a 'finite Weyl curvature' 
singularity, that is to say a singularity at which the Weyl curvature is finite 
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while, of necessity, the Ricci curvature is not. There may be more than 
one plausible way to do this. The strategy adopted here is first to make 
the following definition (after Q): A spacetime {M,gab) is said to admit 
an isotropic singularity if there exists a manifold M D M, a regular Lorentz 
metric gab on M, and a function f2 defined on M, such that 

9ab = ^^9ab for Q> (1) 
n^O on S 

where E is a smooth spacelike hypersurface in M, called the singularity 
surface (note that tilded quantities are defined in the physical space-time, 
while untilded quantities are in the unphysical space-time; also S is smooth 
but there need be no assumption that O is smooth at S). Since the Weyl 
tensor with its indices arranged as 0""^^^^ is conformally invariant and is finite 
in M, it must also be finite in M, so that isotropic singularities form a well- 
defined class of cosmological singularities with finite Weyl tensor. 

Consideration of the conformal Einstein equations near an isotropic sin- 
gularity, for various matter models, leads naturally to a class of singular 
Cauchy problems for the unphysical metric gab with data given on S, for 
which one may seek to prove suitable existence and uniqueness theorems. 
This program was begun in |11] and p!^ where the problem was treated 



in the case of a perfect fluid with polytropic equation of state as matter 
source. There the conformal field equations were written as an evolution 
system for gab and it was shown that solutions which are power-series in 
a time-coordinate could be constructed with data just the 3-metric of the 
singularity surface. The first existence and uniqueness theorem for the class 
of singular Cauchy problems which arise was given by Claudel and New- 
man IQ] for the particular case when the source is a radiation fluid. This 
theorem was adapted to solve the problem for other polytropic fluids in , 
and extended to deal with the massless Einstein- Vlasov equations, (i.e. the 
Einstein equations when the source is massless, collisionless matter satisfy- 
ing the Vlasov or Liouville equation) in Q. There are striking differences 
in the freely-specifiable data in the different cases of perfect fluid and col- 
lisionless matter. For a perfect fluid, the initial 3-metric alone is the free 
data while for the collisionless matter the free data is the initial distribution 
function subject only to a 'vanishing dipole' condition (the initial first and 
second fundamental forms are determined by the initial distribution func- 
tion). It is therefore a natural question to consider other matter models, 
and in particular to see how the gap between these two earlier cases may be 
bridged. 
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In this article, we shall consider a range of other matter models and treat 
them in the style of |11] and |12|. That is, we shall seek to formulate an initial 
value problem with data at the singularity, to identify the data, and to show 
that power-series solutions may be found. We do not seek to manipulate 
the problems into the form in which the Claudel-Newman theorem can be 
applied as this is generally an extremely laborious task. For certainty, this 
would need to be done but in the past the study of power-series solutions 
has been a reliable guide to existence and uniqueness in general. 

The matter models we consider are three: scalar fields, a coupled system 
of Einstein- Yang- Mills- Vlasov studied by Choquet-Bruhat and collaborators 

as a model of a 'quark-gluon plasma', and kinetic theory with either the 
massive Vlasov equation or the massless Boltzmann equations. What we 
find can be summarised as follows: 

For scalar fields subject to a particular condition on the growth of the 
potential with large values of the field it is possible to pose an initial value 
problem for which the data is just the 3-metric of the singularity surface. 
There is no free data for the scalar field, which diverges at the singularity but 
is asymptotically constant in spatial directions, and the second fundamental 
form of the singularity surface is automatically zero. This case is very much 
like the 'stiff' perfect fiuid in|. 

The Einstein- Yang-Mills- Vlasov case is similar to the (massless) Vlasov 
case There are data for the initial distribution function and the initial 
Yang-Mills electric field and magnetic potential. There is a Gauss-law type 
constraint, and a dipole-condition connecting all the data. There are con- 
straints relating the first and second fundamental forms of the singularity 
surface to the matter variables which seem, and in some cases definitely 
are, strong enough to determine them uniquely in terms of the matter vari- 
ables. The massive Vlasov case is essentially the same as the massless case: 
at the singularity, rest-mass of the particles is unimportant. When we in- 
troduce a collision term into the Vlasov (or Liouville) equation to obtain 
the Boltzmann equation however there is a new phenomenon. The collision 
term is defined as usual by a collision integral over products of phase space 
with a kernel determining the differential cross-section. Behaviour of the 
cross-section under conformal rescaling, or equivalently in the limit of large 
energies, now determines the problem. If the cross-section does not grow 
too fast then the initial value problem is like the Vlasov case and the data is 
the initial distribution function, which then determines the initial metric. If 
the cross-section grows fast enough then there is a constraint on the initial 
distribution function, that the collision integral must vanish initially. In this 
case the matter is initially in thermal equilibrium and we are back to the 
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case of a perfect fluid. The data reduces to the initial 3-metric. 

The idea that a study of the Boltzmann equation would lead to an un- 
derstanding of the differences between the perfect-fluid case and the Vlasov 
case is due to Alan Rendall and Keith Anguige [Q, who made some calcu- 
lations in this area, and I am grateful to them for the suggestion. 

The plan of the paper is as follows: we shall give a section to each of 
the three models described above, after ending this section by reviewing our 
conventions for curvature and conformal rescaling. 

Throughout the paper we take the spacetime metric to have signa- 
ture (H ). For a metric gab our definition of the Riemann tensor 

is 

{VaVb-V,Va)V' = R'dabV'' 

while for the Ricci tensor we take 



Rnb — R' 



acb' 



Whenever we consider a rescaling according to (||), tilded quantities will 
always refer to the singular, physical spacetime {M,gab), while un-tilded 
quantities will refer to the regular, unphysical spacetime {M,gab)- For met- 
rics gab and gab related by (|l|), the Ricci tensors are related by: 

Rab = Rab-2VaVbl0gn + 2Val0gnVbl0gn 

-gab (□ log 1^ + 2Ve log logQ) (2) 



and the Weyl tensors by: 



bed — bed' 



If is a unit vector in M then = VLt"" is one in M, and if Wa is a unit 
covector in M, then Wa = Q.~'^Wa is one in M. 

We will usually have a cosmic time coordinate Z at our disposal. This is 
a smooth function in the unphysical space-time which is constant on space- 
like hypersurfaces and vanishes at S. Its gradient is proportional to the unit 
normal Na to the constant-time hypersurfaces: 

Za = VaZ = VNa (3) 

where V'^ = g^^ZaZj,. 

The intrinsic metric hab on the surfaces of constant Z, which is also the 
projection orthogonal to A'"'*, is defined by 

hab = gab - NaNb (4) 
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so that hab is negative-definite. We shall use Da for the intrinsic 3-dimensional 
metric covariant derivative on surfaces of constant Z. 

The second-derivative of Z defines the second fundamental form Kab of 
the constant-time hypersurfaces by the equation 

VaVbZ = V{Kab + NaAb + A^Nb + VzNaNb) (5) 

where Aa = h^^Vb log V = Da log V is the acceleration of the normal congru- 
ence defined by Na, and Vz = V'^N^'VaV (since in general N°-Va = Vd/dZ). 

For the Cauchy problem with data on surfaces of constant Z various 
identities relating 3 and 4-dimensional curvatures are useful. These are 

GabN'^N' = -\{^R-K^ + K^''Kab) (6) 
RacN^K = D,K,^-DbK (7) 

RcdKhl^ = -CNKab + ^Rab + 2KacK^ -KKab 

+DaAb - AaAb (8) 

where Gab,Rab are the 4-dimensional Einstein and Ricci tensors, "^R^^Rab 
are the 3-dimensional Ricci scalar and tensor, K = h'^^Kab is the trace of 
the extrinsic curvature, and ^Cat is the Lie-derivative along A^". When we 
impose the Einstein equations, these identities give the Hamiltonian and 
momentum constraints and the evolution equation for Kab^ respectively. 
The field equations in and were obtained in a first-order form as: 

A\X)^ = A\X)—^ + -B{X)X + C{X, Z). (9) 

Here X is a multi-component vector of unknowns including the first and 
second fundamental forms of the constant-Z surfaces and matter variables, 
Z is the cosmic time coordinate, are comoving spatial coordinates, A^ 
and A^ are symmetric matrices with ^ positive definite, C is analytic in Z 
and B is subject to conditions which will emerge. There is a singularity in 
the time at Z = but, by assumption, nowhere else. 

Suppose we are given data Xq at Z = 0. If we want (^) to hold at 
Z = we at once have a condition, that the term B{X())Xq must vanish. 
We call this the Fuchsian condition for this system. If we now substitute a 
power series in positive powers of Z we can determine the coefficients of the 
series provided none of the eigenvalues of {A^)~^ B{X) is a positive integer. 
With essentially these assumptions, Claudel-Newman [Q] prove existence and 
uniqueness of solution. 
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Putting field equations into precisely this first-order form is tlierefore the 
acid test of existence and uniqueness. However, this is difficult to do. One 
can recognise the Fuchsian conditions and impose the eigenvalue condition 
without going that far and this is what we shall do below. 



2 Scalar fields 

We begin by considering cosmological models whose source is a scalar field. 
We shall find that these may have isotropic singularities provided the po- 
tential satisfies a condition and that there is then an initial value problem 
for which the data is the 3- metric of the singularity surface S. As one might 
expect, this case is similar to the perfect fluid with pressure equal to density 
(see 0). 

If the source of the gravitational field is a scalar field </> with potential 
U{(j)) then the energy- momentum tensor in the physical space-time is 

Tab = Mb - 9ab{^9^'^^c(pd + U{(t>)) 

writing (j)a for Va4>- The conservation equation for this Tab is satisfied by 
virtue of the field equation for the scalar field which is 

a(t> = u'{cf>). (10) 

The tensor which we want to substitute into (^ is the trace-reversed energy- 
momentum tensor: 

Tab - ^fgab = Mb + gabU{4>). (11) 

Now from (^), ( pi] ) and the Einstein equations we find 

Rab = ^VaVbl7-^VaOV60 + 5„fe(^m+^VcJ^V^J^) 

+ ^{^,^, + ^^g^,Um. (12) 

where we don't yet know what to take for $7, but if there is to be an isotropic 
cosmological singularity then this must hold for some and Tab must be 
singular at O = 0. In particular Rab is finite at S, where vanishes, and by 
(|12|) at 0(r2~^) this requires 

~ I (13) 
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in this limit. Looking carefully at (|T2|), we see that we cannot require 17 to 
be smooth at S but we can require Z = Q'^ to be smooth (just as for the 
perfect fluid with p = p). Now to satisfy (|l^) we need </) to be a function of 
Z at least as Z tends to zero and 4> to infinity. Thus near the singularity, (j) 
tends to infinity and the surfaces of constant (p are space-like. Consequently 
we can take a function of (p to be the cosmic time coordinate Z and use this 
to define The correct function to take is determined by (|l^ ) as 

Z = n^ = exp(-) (14) 
a 

where a is a constant to be determined (if a > then (p ~oo as Z ^ 
while if a < then (p — > +c«). We choose units with = 1 and 

substitute (HD into to find 



ZRab = VaVbZ - ^(a' - ^)VaZVbZ + (inZ + Z^U(cP))gab (15) 

where we have also multiplied through by Z. Now we can see that if (H) IS 
to hold at S, where Z = 0, then necessarily = 3/2. Furthermore, with 
this choice we have a regular field equation in the unphysical space-time 
provided the term with the potential U is finite at S. This requires Z'^U{(p) 
to be finite as Z tends to zero, or equivalently exp(^)C/((/>) to be finite as 
(p — > iboo where the correct choice of sign depends on the sign of a. 

If we calculate the contracted Bianchi identity in the unphysical space- 



time using (15) then we find that this is equivalent to the equation 

UZ = -Z^U'U) (16) 
a 



which can also be obtained from (10) by conformally transforming and 
changing the dependent variable. 

We next need to formulate a Cauchy problem for the Einstein equa- 
tions in the form ( p^ ) with the matter equation (|l6|). We perform a 3+1- 
splitting with respect to the hypersurfaces of constant Z and use Z as a 
time-coordinate. The variables in the Cauchy problem will be V, hab and 
Kab at each value of Z and we need evolution equations for each of these, 
with the possibility of constraints between them. 

In terms of comoving coordinates x*, i = 1,2,3 the metric can be written 

as 

dZ"^ 

ds^ = + hijdx'dx^ (17) 
with V and hij as in (^) and (^). 
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From the trace of (|5|) we obtain 

nZ = V{K + ^) (18) 

where K = g°'^Kab, which with (^) gives an evolution equation for V. 
Substituting from (|5|) and (16) into (|l5|) puts the Einstein equations into 
the form 

ZRab = V{Kab + NaAk + AaNh + VzNaNb)+gab{^Z^U'{^) + Z^U{(l))). (19) 

From (^) and (^) we obtain the Hamiltonian constraint as 

- IGabN^N^ = ^R-K^ + KabK"'' = 2(^ + Z[/(0)) (20) 

Zj 

and from (^) the momentum constraint as 

RacN'^K = DcK,' - DbK = ^Db log V (21) 
where we have used the definition of Aa from (^) . For the evolution equations 



we have the evolution of V from (16) and (p^): 

then the evolution of hab, which is essentially the definition of Kab- 

^Nhab = '^Kab (23) 

and finally from (^) and ( p!9[ ) the evolution of Kab 



V 

1 



CNKab = ^Rab + 2KacK{;-KKab + DaAb-AaAb-^Kab 



+hab{—Z^U'{(fy + Z^Um. (24) 
la 



To summarise, the evolution equations are (|2^), ( |2^ ) and (p^), with con- 
straints (pO|) and (^T|). Because the contracted Bianchi identity for (|T5|) is 
( |T6| ) which in turn is ( p2[) and part of the system, we know that the evolution 
preserves the constraints. The evolution equations have a singularity in the 
'time' at Z = 0. In order to put them in the form of (|9|) it would be nec- 
essary to introduce many more variables to make the system first-order in 
space. Rather than do this, we can draw similar conclusions in the present 
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form. Specifically we can read off the Fuchsian constraints and check the 
eigenvalue conditions noted after @. 

For the Fuchsian constraints, if (|2^) is to hold at Z = then we need 
DhV = so that V is constant on S. By a constant rescaling we can suppose 
that V = 1 there. Now for (pO|) at Z = we need initially 



K = - lim{Z^U{(l))) (25) 
while from ( |2^ ) we need initially 

Kab = -habilM^Z^U'icI)) + Z^UicP)). (26) 

For consistency between (25) and ( p6| ) we need to impose 

Z^U{^) + -^Z^U'icf)) ^0 as Z ^0. (27) 

This is in addition to the condition found earlier that Z'^U{(f)) have a finite 
limit as Z ^ 0. Assuming this condition on U{(f)) holds then the Fuchsian 
conditions can be solved to find that the free data at S consist of just the 
initial 3- metric hab- The other variables are V ^ which is one without loss of 
generality, and Kab, which is determined by (p^). 

It is straightforward to construct solutions as power series in Z, if one 
assumes power series in Z for U , and these are unique with the given data 
(this happens because the only singular evolution equation is (^) and the 
singular term has a negative coefficient). While we don't convert the system 
into a form in which the theorem of Claudel-Newman |^] can be applied, 
the existence of unique power series solutions has in the past been a reliable 
guide to the applicability of this theorem. 

It is worth noting that S is necessarily umbilic, in that its second fun- 
damental form is a multiple of its metric, and this multiple is in fact a 
constant. Consequently the magnetic part of the Weyl tensor (see e.g. [^) 
is zero at S and the electric part of the Weyl tensor is proportional to the 
trace- free part of the intrinsic Ricci tensor of E. The initial Weyl tensor is 
finite, as it must be, and is zero if and only if S has vanishing trace-free 
Ricci tensor or equivalently has a homogeneous and isotropic metric. But 
now uniqueness for the Cauchy problem would mean that the cosmological 
model was a Friedman-Robertson- Walker model. Thus with scalar fields we 
have the same situation as with polytropic perfect fluids, namely that if the 
Weyl tensor is zero initially then it is always zero. 
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3 The Einstein- Yang-Mills-Liouville equations 



The Einstein- Yang-Mills-Liouville equations have been considered by Choquet- 
Bruhat and coworkers as a model of the (early) universe where the matter 
content is a 'quark-gluon plasma' [^j. The model has a distribution function 
for massless, collisionless matter with colour-charges which produces a Yang- 
Mills field. With space-time indices a,b,... and Lie algebra indices a,/3, ... 
the matter variables are the Yang-Mills potential and the distribution 
function f{x'^,ph, g") where the vector is a vector of colour-charges. 

We retain the convention that tilded quantities relate to the physical 
space-time and untilded to the rescaled, unphysical space-time. Li fact 

, f and ph are all unchanged under rescaling, as is the Yang-Mills field 
F"^, which is defined by 

F^, = V, A" - V,A," + c"^,^f A" (28) 

where c*^^ are the structure constants for the Lie algebra, say Q. 

The constituent 'particles' of the collisionless matter follow the Yang- 
Mills version of the Lorentz force law which, for a particle with colour- 
charges given by can be written as the system: 



ds 



-he 



= ni.n„-\-n n^. aF'^'" 

ds 

dq°' 



ds 



C /37^a P Q 



here rjap is the metric on Q (assuming it is semi-simple). Again the system 
( p9|) is conformally-covariant in that it takes the same appearance if tildes 
are inserted throughout. 

The Vlasov or Liouville equation for the distribution function follows 
from ( p9D and the chain rule: 

df dx" _^ d£dpa ^ df_d(f_ _ ^ 



dx°- ds dpa ds dq°' ds 

The Yang-Mills current is obtained as an integral over the distribution func- 
tion. Now we do need to watch the conformal rescaling. In the unphysical 
space-time we consider the vector field: 

a = Q Pafu)pUJq 
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where tOp is the standard Lorentz-invariant volume-form on the null cone, 
which can be written with respect to a 3+1-splitting as 



dpidp2dp3 (30) 



and LOq is an invariant volume- form on Q. 
Under conformal rescaling we have 



p 



n-^ujp (31) 



and then the physical Yang- Mills current is 

= n-'j^. (32) 

The physical Yang-Mills equation relates the divergence of the physical 
Yang- Mills field to the physical Yang-Mills current. Now this equation trans- 
forms well under conformal rescaling and in the unphysical variables we find 

VaF"'"' + c^^^Af F'^^T = J*'^ (33) 

The energy- momentum tensor is a sum of two terms. For the Yang-Mills 
field consider 

rpYM _ „ frpa rpPc 1 rpa T?cdf3\ 
-Lab -Val3[^ac^b - ^9ab^ cd^ ) 

and for the collisionless matter consider 



'^ab = j fPaPb^^p^^q- 



Note that both these tensors are trace-free. The physical energy-momentum 
tensor is then 

f,, = n-\T^,^'' + T^,). (34) 

We write down the Einstein equations in the unphysical space-time with 
the aid of (^. Following the argument in [|| we may choose the cosmic 
time coordinate Z so that n\Z = everywhere and y = 1 at S and we can 
then take Z = Vt (there is sufficient conformal gauge freedom to allow this, 
provided there is an isotropic cosmological singularity, because the energy- 
momentum tensor is trace- free). From (|2|) we obtain 

Z^Rai, = 2ZVa^bZ - AVaZVbZ + V^gab + T^b^' + T^b- (35) 

The dependent variables on each surface of constant Z are V, hab and Kab 
as in section 2, together with the Yang-Mills potential and the Yang- 
Mills electric field := F^^N^ , and the distribution function /. There 
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is sufficient freedom in the conformal gauge to ensure that K = at T,. We 
choose the temporal Yang-Mills gauge so that N°'A" = 0. The Yang- Mills 
equations ( p3| ) imply a Gauss-law constraint: 

which translates as 

D^E'^" + c%Al^E''^ = j q"fdpidp2dp3. (36) 

This constraint holds on every surface of constant Z and in particular needs 
to hold on the data at S. 

We perform a 3+1-splitting with respect to the hypersurfaces of constant 
Z as before, and decompose the Einstein equations ( |35| ) into evolution and 
constraints. The Fuchsian constraints for data at S will follow from the 
condition that the right-hand-side in (35) must be 0{Z'^). We write spatial 
indices i,j,k,... = 1, 2, 3 and then the conditions arise from the vanishing 
of the 0(1) terms in the and (0, i) components in (|35|) and the van- 
ishing of the 0{Z) terms in the components. (The 0{Z) terms in the 
(0, i) components vanish by virtue of these conditions since the contracted 
Bianchi conditions hold.) The conditions are most easily written in terms of 
multipole moments Xi, Xijj Xijk, ■ ■ ■ of the matter density and J", J^", J^"^, . . . 
of the charge density which we define at S by 

1 

Xi = 
Xij = 

Xijk = / , Q 




and so on, and 



1 /" <f!Vi ,3„ 



. r^h J (p0^2 



ja _ _J_ [ (I'^fViVlVk ,. 



_ 1 /" <ffPiPjPk^,3 
'^^^ ~ (p0)3 

and so on, where it is understood that f,h and are evaluated at S. We 
have used (|30|) and in these definitions 
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Now the Fuchsian conditions are 



= h,,+x^J-vMErE^ -KkFp 

+\h,{2EtE'^P + F^^F'^P)) (37) 
= X'.-riapE^'^F^^ (38) 



and the most complicated which can be written 



-r,.,(4i?^ii;^ + 2F(7"jj5„ + i^^"4 

-DkE^F^^^ - D^,E^k\F;f + 2EID^F^^^ 

+ hi,{b\Fkrr,^Ei) - E^^4)) (39) 



where 



Qlmij — ^^lif^mj hlmhij Xlmij 

-vME^^E^hmj + E^E^hrm + E^E^^hi, + E^E^hu 
—EfE^him — EfE^hij 

-Ff^F^ - Ff^ f'^ -i- h- Fff F''f^ -i- h, F^ F^^^ 

^il^jm ^ im^ jl ^ "'■I] ^ Ik^ m ^ 'Hm^ ik^ j 

-huhm,{2EtE^P + F^^F^^^) 

+himh,{\EtE^P - \Fi^F^-^)). (40) 



If we consider first (|3j) , then we recall from |2| that in the case of vanishing 
Yang-Mills fields this actually determines the initial metric hij from the 
initial distribution function /. It is hard to see whether this remains true 
for arbitrary Yang-Mills fields but it will certainly remain true for small 
fields. 

Next note that, by virtue of (|37|) and (|38|), Qimij is symmetric and trace- 
free on each index pair, and symmetric under interchange of the two pairs. 
If we think of it as a map from trace-free symmetric tensors to trace-free 
symmetric tensors, it necessarily has real eigenvalues and we can solve (39) 



for Kij provided none of these eigenvalues is zero (recall that Kij is trace- 
free at E by choice of conformal gauge). Again this is a difficult question in 
general although we know from [|[ that it is true if the Yang-Mills fields are 
zero, so that it will be true if they are small. 

Finally ( pSj ) and ( ^6] ) are constraints on the initial matter variables. After 
the Fuchsian conditions we should check the eigenvalues as in the discussion 
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following (|9|) , but again we can argue that they must still satisfy the required 
condition, at least for small fields. 

To summarise, the data at S consists of the Yang-Mills potential Af, 
from which the magnetic field F^^ is calculated by (^), the Yang- Mills elec- 
tric field , the distribution function / and the first and second fundamen- 
tal forms hij and Kij. These are subject to the Gauss-law constraint (]36| ) 
and the three Fuchsian constraints (|37|), (|38| ) and (|39|). In the pure Vlasov 
case, with the Yang-Mills fields omitted, (^) suffices to determine hij from 
/ and (|39| ) then determines K^j from / (see or for details of this). 
Thus in that case the free data is just / subject to the counterpart of (p^. 
With the Yang-Mills fields present it is not so easy to see that (l37|) deter- 
mines a unique hij for arbitrary Yang-Mills fields and that (39) determines 
a unique Kij, but this must certainly be true for small fields by continuity. 
Once the Fuchsian constraints are satisfied, and relying on the experience 
of earlier cases, we expect to have existence and uniqueness of solutions. 



4 Massive Einstein- Vlasov and Einstein-Boltzmann 

In this section we do two things. First we look at massive Einstein- Vlasov in 
the expectation that near the singularity there should be no difference be- 
tween massless collisionless matter and massive collisionless matter. Then, 
restricting to massless for simplicity we seek to drop the collisionless condi- 
tion. This is partly to move to a more realistic situation but also to resolve 
a puzzle left over from earlier work and noted in the Introduction. 

For an isotropic singularity with a perfect fiuid (at least for linear equa- 
tions of state 1^) the data are just the initial 3-metric. There are no separate 
data for the matter. For an isotropic singularity with massless collisionless 
matter the data is just the initial distribution function, which determines 
the initial first and second fundamental forms This time there are no 
separate data for the geometry. It was suggested to me by Alan Rendall [^] 
that by including a collision term in the Vlasov equation, in other words 
by turning to the Einstein-Boltzmann equations, one might find a bridge 
between these two extreme examples. This turns out to be the case. 

We begin by considering the massive Einstein- Vlasov equations, by which 
we mean the Einstein equations whose source is collisionless matter where 
the matter is now supposed to be composed of particles of a single nonzero 
mass, say m. The presence of the mass changes the detail of the equations 
but we shall see that that change has no infiuence near to the singularity. 

In this case, the stress-energy-momentum tensor in the physical variables 
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IS 



rpL 



J fPaPbOJp (41) 



with ^ 

= ,Q ^^ dpidp2dp3. (42) 



PV-S' 

The distribution function / is supported on the mass shell where 

g^^PaPb = 



which we must solve to find p^ in (^2|) . We subsequently regard p^ (and later 
p*^) as a function of pi and the coordinates, and the distribution function as 
a function only of pi and x"" . 

Under conformal rescaling according to (||) we still have (|3^) with ojp as 
in (^) but now with 

pO = v^po = V{-h'^piPj + m^rj^)! . (43) 

As one might expect, the mass becomes insignificant near the singularity. 



where Q vanishes. The unphysical energy-momentum tensor to go into (35) 



may be written, with the aid of the metric (17) and the expression (H3) as 



L 



fPaPbdpidp2dp3 



V^J {-h^'^pkPm + m?^"^)^ 

Provided / is supported away from the origin in p-space, which one usually 
assumes for regularity (see e.g. |l|]), the energy momentum tensor can be 
expanded as a series in positive powers of In particular in the limit at S 
it is the same as in the massless case. Consequently in rescaling we take Vt 
to be Z just as in the massless case. 

We also have to change the Vlasov equation in the presence of mass. The 
Vlasov equation expresses the vanishing of the derivative of the distribution 
function along the geodesic flow, or equivalently the vanishing of its Poisson 
bracket with the Hamiltonian. In the unphysical variables this is 

C^f[Pa.x'') = {\r'paPbJ}pB 

— 9 Pa r, I, no^- PaPb 

ox" 1 OX'^ Opc 

= 0. (44) 
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Written in this form, it is easy to see how to conformally transform the 
equation. With the metric as in (|l^ and rescahng as in (|l]) we obtain 

The last term in this actually vanishes because is Z, but there are still 
explicit appearances of Z in ( ^5|) through p^. However these are always in 
the numerator and we get no new singularities from them. 

We conclude that, at least at the level of looking for power series solutions 
in the time- variable, there should be no difference between massive and 
massless Vlasov. 

To go from Vlasov to Boltzmann is to go from collisionless matter to the 
inclusion of a collision-term in (^). Symbolically we have 

C^f = C{fJ) (46) 

where C(/, /) is a bilinear functional in /. We set up the formalism for this 
following Choquet-Bruhat We restrict to 2-body collisions and suppose 
that incoming particles with 4-momenta p^^^ and p^^^ (temporarily suppress- 
ing space-time indices) collide to give outgoing particles with 4-momenta 
p^"^^ and p^^\ where all particles have mass m and total 4-momentum P is 
conserved: 

p = p(0)+p(l)=p(2)+p(3) 

If we choose p^^^ then there are 3 degrees of freedom in p^^-*, two in p^"^^ and 
none in p^^^ . The collision term is written as an integral: 

C{f,f) = I (/2/3 - /o/i)A;(^",P^°\p^'^P^'\p^'^H(i)6 (47) 

where fn = f{x°',Pi^^), A; is a kernel to be discussed below, ojp is as in ( pO| ) 
but with the time-direction defined by p^^^ , and .^2 is an invariant (Leray) 
2-form on the space of allowed p^^^ (see [^] for a precise definition). The 
kernel k is related to the differential cross-section and is usually assumed, 
for reasons of symmetry, to be a function only of the total energy involved 
in the collision, which is 

s = r'PaPb (48) 

and the angle 6 between p^^^ —p^^^ and p^"^^ —p^^^ , which is the scattering angle 
in the centre-of-mass frame, so we may write k{s,9). For Choquet-Bruhat 
A; is the differential cross-section but for other authors, e.g. Grout et 
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al Q, k = sa{s,9) where a is the differential cross-section. The important 
consideration for us is the behaviour of the various terms in ( ^7|) under 
conformal rescahng. We know how to transform s from (^ ) and Op from 
(|3l|). The Leray form ^2 is conformaUy-invariant ^ and so the rescahng of 
(^6|), at least for the massless case is 

Cxf = ^''C^f 

= J (/2/3 - fofi)m-^s, ep^^a- (49) 

We want the Boltzmann equation ( p9| ) to hold up to S and so we need to see 
if it gives another Fuchsian condition. This will depend on the behaviour of 
the collision term and specifically on the behaviour of the kernel k as Q goes 
to zero or equivalently on the behaviour of k(s, 6) as s increases without 
bound. 

If k is sufficiently well-behaved that the integral in ( ^9[ ) exists in the limit 
for every / then we get no new Fuchsian condition. The situation is just 
like the Vlasov case dealt with in ||^] and we should expect existence and 
uniqueness with /q as the free data (subject to the Fuchsian constraint cor- 
responding to ([3^), that Xi = 0)- The other Fuchsian constraints determine 
the first and second fundamental forms of S in terms of /q. 

However if k diverges for large s then we do get a new Fuchsian condition: 
for the collision term to be finite, the term (72/3 — /o/i) in the integral must 
vanish at S. As is familiar from the literature (see e.g. |^]), this condition 
requires the initial distribution function to take the form 

h{x\pj) = exp(-a(x*) - p''{x')pa) (50) 

for a function a(x*) and a vector- field /3"(x*), both defined at S. Next, the 
Fuchsian condition = when imposed on a function /q of the form of (|50| ) 
requires that the vector be proportional to the normal A^" to S. Now 
the initial distribution function takes the form 

/o = exp(-a(x^) - /3(x^)po) (51) 

which is formally the distribution for local equilibrium with a and (3 related 
to chemical potential and inverse temperature (for global equilibrium these 
would be constants). Note that j3 is not the physical inverse temperature. 
This would be /3 = Vtf3 with the rescahng conventions in use here, and this 
vanishes at S as one would expect. 
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The quantity po in ( [5l| ) is understood as a function of and pj as in 
(H) but now 

since at S we can take V = 1 (and we are supposing for simplicity that 
m = 0). With this in (^), the Fuchsian constraint corresponding to (|3^ ) is 
satisfied provided a and /3 are related by 

^4g-« = 2567r2. 

Finally the third Fuchsian constraint, corresponding to (|39|), forces the initial 
second fundamental form to vanish, because the third moment Xijk of /o 
vanishes for a distribution function of the form of (|5l|). Thus the data 
consist of the initial metric hij and a function which we can call the 

(unphysical) initial inverse temperature. The matter is initially in local 
equilibrium, so that the stress-tensor is the stress-tensor for a perfect fluid. 
However, there is no reason to suppose that this form for / will be preserved 
by the evolution 



In conclusion, we see that, as suggested by Rendall and Anguige the 
Boltzmann case sits between the perfect-fluid and Vlasov cases studied in 
and Q, and depending on the behaviour of the differential cross-section at 
large energies it will, as regards data at the singularity, approach the Vlasov 
case or the perfect fluid case but with an extra function. 
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